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O ■ ABSTRACT 

' Using far (FUV) and near (NUV) ultraviolet photometry from guest investigator pro- 

\ grammes on the Galaxy Evolution Explorer (GALEX) satellite, optical photometry 

_ • from the MDM Observatory and optical integral-field spectroscopy from SAURDN, we ex- 

t— | ' plore the UV-linestrength relations of the 48 nearby early-type galaxies in the SAURDN 

sample. Identical apertures are used for all quantities, avoiding aperture mismatch. 
We show that galaxies with purely old stellar populations show well-defined correla- 
tions of the integrated FUV— V and FUV— NUV colours with the integrated Mg6 and 
5_J , H/3 absorption linestrength indices, strongest for FUV— NUV. Correlations with the 

NUV— V colour, Fe5015 index and stellar velocity dispersion a are much weaker. These 
correlations put stringent constraints on the origin of the UV-upturn phenomenon in 
early-type galaxies, and highlight its dependence on age and metallicity. In particular, 
despite recent debate, we recover the neg ative correlation betwe en FUV— V colour and 
Mg linestrength originally publicised by iBurstein et all Jl988), which we refer to as 
the "Burstein relation" , suggesting a positive dependence of the UV-upturn on metal- 
licity. We argue that the scatter in the correlations is real, and present mild evidence 
that a strong UV excess is preferentially present in slow-rotating galaxies. We also 
demonstrate that most outliers in the correlations are galaxies with current or recent 
star formation, some at very low levels. We believe that this sensitivity to weak star 
formation, afforded by the deep and varied data available for the SAURDN sample, ex- 
plains why our results are occasionally at odds with other recent but shallower surveys. 
This is supported by the analysis of a large, carefully-crafted sample of more distant 
early- type galaxies from the Sloan Digital Sky Survey (SDSS), more easily comparable 
with current and future large surveys. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
photometry - galaxies: stellar content - galaxies: structure - ultraviolet: galaxies 
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1 INTRODUCTION 



Far-ultraviolet (FUV) radiation was first discovered in early- 
type galaxies by the Or biting Astronomical Observatory-2 in 
1969 (ICode et al.lll972T ). It was rapidly realised that, below 
« 2000 A, many early-types show sp atially extended ris- 
ing flux for decreasi ng wavelengths fe.g. lCode fc Welchll979l : 
iBertola et aI]|l98Ch , a surprising behaviour for presumably 
old stellar populations. This spawned a debate on the origin 
of this UV upturn that continues to this day. 

The leading hypothesis suggests that this flux is 
due to hot stars on an extreme horizontal branch 
(EHB) and their progeny (e.g. asymptotic giant 
branch (AGB) stars manque), either metal-p oor (e.g., 



Lee. Demarque fc Zinnl Il994fc [Park fc Led 1 19971; see also 
Buzzoni fc Gonzalez-Lopezliral 20081) or metal-rich (e.g. 
Greggio fc RenzmT 1990l; Horch. Demarque fc Pinsonneauftl 



Bressan. Chiosi fc Fagottc 1994 



1992; 



Dorman. O'Connell fc Rood! Il995l ; lYi. Demarque fc Kiml 



LiK){\ ) , witn possible s ignificant crtccts irqm_tle ennancement 
(e.g. iLee et alj l2005bl ; iKavirai etail 120071 ). Binary star 
systems offer another compelling poss ibility, with limited 
dependence on age and metallicity (e.g. Maxted et alj 200 ll ; 



lHan et all l2003l ; 1rlan. Podsiadlowski fc Lvnas-Gravl |2007j ). 



Post-AGB (PAGB) stars in the planetary nebula stage 
also emit FUV radiation, but they are generally thought 
to be too hot, too few and to have too short a lifetime to 
explain the observed fluxes alone (e.g. iGreggio fc Renzinil 
Il990l : iDorman et ai1ll995l ; iBrown et al.lll997i ). Clearly, very 
young stars in star-forming galaxies can also emit i n the 
FUV (see, e.g., iFerreras fc Silkl |2000| ; lYi et alJliobH ), but 
they necessarily have evolving properties (e.g. characteristic 
temperature of the UV spectrum), they have different spec- 
tral signatures from those normally observed in UV-upturn 
galaxies (e.g. strong C IV lines) and they are expected to 
be clumpy, none of which is supported by obser vations. A 
thorough review of these issues is available in lO'Connelll 
(I1999T ). while iBrownl (|2004l ) an d IyH (120081) p rovide more 
recent but shorter overviews. See iHeberl (|2009l ) for a review 
of hot subdwarf stars. 

In add ition to its roughly constant spectral slope (T e s « 
20,000 K; lBrown etall [19971. one of the main constraints 
on the possible origins of the UV upturn is its systematic 
de pendence on b ulk galaxy properties . Although mentioned 
in (Fabed i|l983l) , iBurstein et aLT i| 19881 . hereafter IBBBFL88I ) 
are generally considered the first to have systematically in- 
vestigated the dependence of the UV upturn on the stellar 
population and dyn amical pro perties of a large sample of 
early-type galaxies. [BBBFL88I found clear non-linear cor- 
relations of the UV-optical colo urs with both the absorp- 
tion linestrength index Mg2 fsee lFaber et al.lll985l ) and the 
central stellar velocity dispersion a. In memor y of David 
Burstein and in view of the long-lasting impact of BBBFL88, 
we will hereafter refer to the former correlation as the 
"Burstein relation" . These correlations have formed the ba- 
sis of much of the observational and theoretical work on the 
UV upturn to date. Crucially, contrary to broadband optical 
and near-infrared colours, the Mg2 linestrength was found 
to be stronger for bluer UV-optical colours, i.e. galaxies be- 
come bluer with increasing metallicity in the UV. Globular 
clusters show different correlations, poi nting to a slightly dif- 
ferent origin for the UV radiation (e.g. IDorman et al.lll995l ; 
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Sohn et al.ll2006l ; but see also iLee. Gim fc Casetti-Dinescul 
20071 ). 



The recent availability of large-area, deep UV imaging 
from the Galaxy Evolution Explorer (GALEX) satellite has 
added another layer of complexity to the problem. Observa- 
tions of large samples of early-type galaxies from the Sloan 
Digital Sky Survey (SDSS) have revealed that a very large 
fraction of them (over 30%) have rem arkably s t rong near- 
ultraviolet (NUV) and FUV emission. lYi et alj l|2005l ) ini- 
tially concluded that the majority of UV-blue early-type 
galaxies were inconsistent with exhibiting a classic UV up- 
turn, and more likely expe rienced residual star formation. 
ISchawinski et all (|2007h and lKavirai et all (120071 ) elaborated 
on this with larger samples, the main conclusion being that 
residual star formation of a few percent of the stellar mass is 
very common among early-type galaxies. The UV emission 
regularly found in early-type galaxies is simply too powerful 
to be accounted for by any current model of the UV upturn, 
and is far more powerful than anything found in traditional 
UV-upturn galaxies such as NGC 1399 and NGC 4552. 

Some recent work has also thrown the results of 
BBBFL88 into doubt. In particular, for a very large sample 
of local early-type galaxies with GALEX U V imaging and 
SDSS optical spectroscopy, iRich et all l|2005h report no cor- 
relation of the FUV— r colour against either the metallicity- 
sensitive Mg2 and D4000 indices or the stellar velocity dis- 
persion, even after attempting to remove non-early-type and 
star-forming galaxies from their sample. If anything, a slight 
correlation opposite to the Burstein relation is observed. 
iDeharveng. Boselli fc Donas! l|2002l ) report similar findings, 
although without active galactic nucleus (AGN) or star- 
formation cuts and using a UV filter centred on 2000 A, 
where the UV upturn phenomenon may not dominate. Us- 
ing data on nea rby galaxies but so rting according to mor- 
phological tvpe. lPonas et all |2007h do find correlations for 
elliptical galaxies but those for lenti culars are weak at b est, 
with much scatter and many outliers. iBoselli et al.l (120051) re- 



BBBFL8S for 



cover weak correlations analogous to those of . 
Virgo Cluster early-type members, but it is unclear whether 
the correlations extend to low-luminosity systems. If these 
recent results prove to be robust, much of the theoretical 
work on the UV upturn will need to be revised, particularly 
its dependence on metallicity. 

To sort out possible discrepancies between the above, 
as much information as possible is required on the sam- 
ple galaxies under consideration, particularly regarding pos- 
sible low-level star formation. We thus consider two sam- 
ples here. First, the SAUR0N sample of 48 nearby ellipti- 
cal a nd lenticul ar galaxie s (see Table [T] Ide Zeeuw et al.l 
2002, hereafter IPaper III) , for whi ch detailed integral- 
field map s of the st ellar kinematics fEmscllcm et al 1 l2004l . 
hereafter IPaper III), stellar absorpti on linestre ngths and 
populations dKuntschner et al.l |2006|. hereaft er IPaper Vll ; 
iKuntschner et al.l 120101 ." hereafter IPaper XVlJ) and ionised 
gas d istributio n and kinematics ( Sarzi et all |2006l. here - 
after IPaper Vt ISarzi et al.l l20ld . hereafter Paper XVj) 
available, along with d etaile d dynamical modeling 



arc available, along_^ witn d etaile d dynamical modenng 
(e.g. ICappellari et all l2006l, 12007); lEmsellem et all 120071 , 



hereafter Paper ivl ~ Paper Xl and Paper IX, respectively) 



The molecular gas conte nt (| Combes. Young fc Bureau 



n r 



2007; ICrocker et all 120101 ). UV emission (|Jeong et al 
2OO9J, hereafter IPaper XIlJ) and mid-infrared (MIR) emis- 



sion dTemi, Brighenti fc Mathewsll2009l ; IShapiro et aill2010l . 
hereafter Paper Xvl) of the SAUROH sample have also been in- 
vestigated, to more directly probe its star formation proper- 
ties. Secon d , we c onstr uct a SDSS s a mple similar to those of 
IRich et all l|2005l ) and iDonas et all |2007l ). but apply much 
more stringent cuts to remove star-forming galaxies and low 
quality data. In bot h cases, we recover the relationships orig- 
inally suggested bv lBBBFL88l for the non-star-forming (i.e. 
quiescent) galaxies only, while including galaxies with even 
a very low level of star formation considerably blurs the re- 
lations or can even give rise to trends in the opposite sense. 

In Section [2] we present the photometric and spectro- 
scopic data and related analysis. Section[3]presents the UV- 
linestrength relations for the SAUR0N sample while Section U 
discusses the effects of star formation on them. The SDSS 
sample is presented and discussed in Section [5] We bring 
all our results together in Section HH and conclude briefly in 
Section [7] 



2 OBSERVATIONS AND DATA REDUCTION 

Although spatially-resolved data provide in principle a more 
powerful discriminant of UV-upturn theories, we focus here 
on integrated qu antities as they provide a link with past 
work (particularly BBBF188) and offer a crucial benchmark 
for future studies at high redshifts. The latter in particular 
are essential to probe the evolution of the UV-upturn with 
look-back time (and thus turnoff mass), an i mportant test of 
single versus binary stellar origin s (see, e.g.. lYi et al.lll999l ; 
ILee et al.ll2005al : iRee et all I2007T ) . The two-dimensional na- 
ture of the SAUR0N data will be more fully exploited in a 
com panion pa per (Jeong et al., in preparation). 

IPaper "vl showed that, to first order, iso- index contours 
follow a galaxy light (although they are slightly flatter than 
the isophotes for Mg6 in ~ 40% of the cases). We therefore 
adopt elliptical apertures to measure integrated quantities 
throughout this paper. We will generally refer to quantities 
(magnitudes and linestrengths) integrated within ellipses of 
semi-major axis Re/2, where R e is the effective radius of the 
galaxy (i.e. the radius of the aperture encompassing half the 
galaxy light). A radius of R c /2 is ideal for our work as it 
is large enough to accurately represent global values, yet it 
is small enough that the corresponding elliptical apertures 
are fully covered by the SAURON field-of-view (FOV) for most 
objects (~ 80%). Other interesting radii include R e /8 (i.e. 
central values) and R e (global values), although the latter 
requires a small extrapolation for ~ 60% of the objects. 



2.1 UV photometry 

All 48 SAURON early-type galaxies except seven were observed 
in FUV (1350-1750 A) and/or NUV (17 50-2750 A) with the 
Medium imagin g mode of GALEX Csee lMartin et al.ll2005l ; 
iMorrissev et al.l 120051 ) . A total of 17 galaxies were observed 
as part of our own UV imaging survey of the SAURON sam- 
ple (GALEX guest investigator programmes GI1_109 and 
GI3_041) and another 17 as part of the GALE X team. Nearby 
Galaxy Survey (NGS; lGil de Paz et al l [20051 ') . All these ob- 
jects were observed for at least one orbit (typically 1500- 
1700 s) and have excellent data for any aperture up to R c . 
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Table 1. SAURQN sample data 



Galaxy FUV flc/2 NUV flc/2 V Rc/2 E(B-V) Mgb Rc/2 Fe5015 flc/2 IL3 He , 2 a c Survey 

(mag) (mag) (mag) (mag) (A) (A) (A) (km s — 1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
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Columns: (2)-(4) and (6)-(8): Quantities integrated within elliptical apertures of semi-major radii R c /2; (2)-(4) have been corrected for 
Galactic extinction using (5) and the errors are based on Poisson statistics and our sky subtraction. Th e adopted for mal random error 
on the values in (6)— (8) is 0.1 A (see lPaper VTTl , (9) Effective stellar velocity dispersion within Re from I Paper XVlH . with an adopted 
error of 5%. (10) Origin of data: A=All-sky Imaging Survey, G=own guest investigator programmes, N=Nearby Galaxy Survey. 



The remaining 7 galaxies have only short exposures (typ- 
ically 100-150 s) from the GALEX team All-sky Imaging 
Survey (AIS; see lMartin et al]|2005l ). and magnitudes inte- 
grated within R a are less reliable. The resulting sample is 
listed in Table [1] The data reduction, analysis and product s 
are described more fully in a separate paper (|Paper Xllj) . 
so we only summarise our processing here. 

We used the pre-processed images delivered from the 
GALEX pipeline but performed our own sky subtraction by 
measuring the background level in source-free regions of the 
images. The spatial resolution is roughly 4'.' 5 and 6'.'0 FWHM 
at FUV and NUV, respectively, sampled with l'.'5 x l"5 pix- 
els. To avoid spurious colour gradients (especially in the in- 



ner parts), we thus convolved the FUV data to the spatial 
resolution of the NUV observations before any analysis. 

We carried out surface photometry by measuring the 
surface brightness within elliptical annuli in the standard 
manner, using the ELLIPSE task in IRAF (Image Reduc- 
tion and Analysis Facility). To derive accurate colours, the 
ellipses were fitted to the NUV images only (with the bet- 
ter signal-to-noise ratio S/N), and those ellipses were then 
imposed on the FUV image s. Pho tometric zero-points were 
taken from iMorrissev et alj (|2005) and the data were cor- 
rected for Galactic extinction usi ng the extinction law of 
ICardelli, Clayton fc Mathlsl (jl989t ). We take the E(B - V) 
value for each galaxy as tabulated in the NASA/IPAC Extra- 
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galact ic Database (NED) from lSchlegel. Finkbeiner fc Davisl 
(|l998h - For this paper, only the FUV and NUV surface 
brightness profiles integrated within R a /2 are discussed, but 
we have checked that analogous results stand with apertures 
of Rq/8 and R c . The integrated UV magnitudes corrected 
for Galactic extinction are listed in Table [TJ along with er- 
rors based on Poisson statistics and our sky subtraction (all 
UV magnitudes are AB magnitudes). 



2.2 Optical photometry 

Optical images in the Hubble Space Telescope (HST) fil- 
ter F555W (similar to Johnson V) were obtained for the 
entire SAURON sample using the MDM Observatory 1.3-m 
McGraw-Hill Telescope over several observing runs in 2002- 
2005. The field-of-view of the MDM images is 17.'3 x 17'.3 
with / .'508 x 0'.'508 pixels, permitting excellent sky subtrac- 
tion and proper sampling of the seeing. The MDM observa- 
tions are described in Falcon- Barroso et al. (in preparation) . 
They were reduced in the standard manner in IRAF. The 
photometric calibration wa s obtained by observing photo- 
metric standard stars from lLandoltl (|l992h over a range of 
air masses, resulting in a systematic calibration uncertainty 
of 0.03 mag. Like the FUV data, the MDM images were con- 
volved to the spatial resolution of the NUV data before any 
analysis. 

Again like the FUV data, the surface brightness profiles 
for the MDM data were measured only after imposing the 
NUV ellipses, ensuring accurate UV-optical colours. Mag- 
nitudes integrated within elliptical apertures of semi-major 
axis R a /2 (identical to the UV apertures) were then calcu- 
lated and corrected for Galactic extinction. The integrated 
and corrected V magnitudes (and errors) are listed in Ta- 
ble [U 



2.3 SAURON linestrengths and kinematics 

All linestren gth measu rements co me from the S AURON data 
published in lPaper Vll updated in IPaper XVI j The obser- 
vations and data reduction were thoro ughly discu ssed there 
and will not be repeated here. While 'BBBFL88 strived to 
derive integrated V magnitudes within apertures equivalent 
to the roughly 11" x 22" rectangular aperture of the Interna- 
tional Ultraviolet Observer (IUE), they had no choice but to 
use essentially central measurements from slit spectroscopy 
for Mg2 and a. A significant improvement in the current 
work is that integral-field spectroscopy allows us to measure 
linestrength indices over exactly the same apertures as the 
UV and optical p hotometry. 

IPaper XvTj] lists linestrength values integrated within 
circular apertures of radii R c /8 and R c , but we again use 
here NUV-defined elliptical apertures of semi-major axes 
Re/2, as done in Sections l2.ll andl 2.2l to avoid potential aper- 
ture mismatch problems. We use simple luminosity-weighted 
linestrength averages of all the bins within the apertures. 
The resulting integrated linestrength values are listed in Ta- 
ble]]] While not formally equivalent, we show in Appendix IA1 
that these measurements are as robust as those obtained by 
summing the spectra of all the bins concerned and rederiving 
the linestrengths from scratch. We therefore ad opt the same 
uncertainty for our integrated linestrengths as lPaper XVlj . 



a formal random error of 0.1 A (maximum) for all measure- 
ments. 

Since linestrengths are known to cor relate with the 
central velocity dispersion of galaxies (e.g . Terlevich et al.l 



19811; iBurstein et al.lll988l ; iBernardi et alTl2003l ), following 



BBBFL88 we also explore the dependence of the UV 



optical colours on a. For that, we do not use central veloc- 
ity dispersion measurements as is usually done, but rather 
again use properly integrated measurements obtained by 
summing all the spectra within a given aperture, and re- 
measuring the stellar kinematics (assuming a pure Gaussian 
line-of-sight velocity distribution (LOSVD), i.e. no high or- 
der Gauss-Hermite terms). Such stel lar velocity dispersion 
measuremen ts are already available in lPaper IVl (updated in 
IPaper XVIl|) for a circular aperture of radius i? c , and the ex- 
trapolation required when the field-of-view of SAURON does 
not fully cover R a is detailed there. The resulting "effec- 
tive" stellar velocity dispersion <r is truly the luminosity- 
weighted second moment of the LOSVD within _R C , and is 
a good approximation to the second velocity moment which 
appears in the virial equation. Integrated velocity disper- 
sions only weakly depend on the aperture used, so we do 
not re-calculate new values for the elliptical apertures used 
throughout this paper, bu t instead sim ply adopt the values 
and uncertainties listed in IPaper XVlj (see Table [TJ). 



3 UV RELATIONS 

When exploring possible correlations between stellar ab- 
sorption linestrengths and broadband photometr y, our main 
advantage over all p revious studies, inc luding BBBFL88, 
iRich et al.l (120051 ) and lDonas et all l|2007l '). is the availability 
not only of UV imaging, but most importantly of integral- 
field spectroscopy for all targets. All photometric and spec- 
troscopic measurements are thus made on identical aper- 
tures, as described in Section [2] 

Figure [TJ shows the distance-independent correlations 
between the FUV— V colour and the linestrength indices 
Mgfe, Fe5015 and H/3. Although the transformation of Mgb, 
Fe5015 and H/3 to stellar age, metallicity and Q-element 
(over-) abundance is model dependent, H/3 is primarily a 
tracer of age. For the range of metallicities present in early- 
type galaxies and c onsidering various st e llar population syn- 
thesis models (e.g. I Thomas et al.l [20031 ; ISchiavonll2007t ). an 
H/3 linestrength of 1.8 A corresponds to single stellar popu- 
lation luminosity-weight ed ag es of roughly 5 to 10 Gyrs (see, 
e.g., Figure 3 in IPaper XVlj) . This linestrength value thus 
represents a good choice to separate young from old galax- 
ies, or rather primarily old galaxies with a significant amount 
of recent star formation and young stars from primarily old 
galaxies with no or a minimal additional young stellar pop- 
ulation. Changing this threshold value slightly does not sig- 
nificantly affect any of our results. To highlight the effects of 
recent star formation in Figure [TJ galaxies with H/3 > 1.8 A, 
that are inconsistent with a purely old stellar population, are 
therefore plotted as crosses. Galaxies with no or minimal 
(in a luminosity-weighted sense) young stellar populations 
(H/3^ 1.8 A) are plotted as filled circles. 

There is no physical reason why a linear relation should 
exist between linestrengths and UV excess, but linear fits to 
the data with H/3^ 1.8 A are shown as solid lines to guide 
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Table 2. Parameters of the best-fit linear UV-linestrength rela- 
tions for old galaxies only 



Colour 


Line 


Slope 


Zero-point 


Scatter 


(mag) 


(A) 


(A mag -1 ) 


(A) 


(A) 


FUV-V 


Mgfe 


-0.40 ±0.09 


6.75 ±0.57 


0.15 




Fe 


-0.07 ±0.11 


5.11 ±0.75 


0.22 




H/3 


0.12 ± 0.05 


0.83 ±0.32 


0.08 


FUV-NUV 


Mgfe 


-0.64 ±0.12 


4.89 ±0.14 


0.15 




Fc 


-0.21 ± 0.17 


4.89 ±0.19 


0.21 




H/3 


0.20 ±0.08 


1.40 ±0.09 


0.07 



the eye. The value of the correlation coefficient is reported 
in the bottom right corner of each panel, and the parame- 
ters of the best-fit lines are reported in Table [2l all taking 
errors properly into account. Lastly, the symbols are colour- 
coded according to the effective stellar veloci ty dispersion 
cr c , increasing from purple to red (|Paper XVI J) . 



3.1 Burstein relation: FUV-V vs. Mgfe 

The top panel of Figure Q] shows the (FUV-V)-Mg b rela- 
tion, analo gous to the (1550— V)-Mgz relation discussed by 
IBBBFL88I . although both the FUV filter and the Mg in- 
dex are slightly different and we use a common aperture of 
Re/2. This is thus our own version of the Burstein relation. 
Reassuringly, we see a very similar trend. Most importantly, 
there is a clear correlation between the FUV— V colour and 
the Mg& linestrength for galaxies with no significant young 
stellar population (H/3 ^ 1.8 A, i.e. quiescent galaxies). This 
is the systematic variation of t he hot comp onent of old stel- 
lar populations discovered by IBBBFL88T Galaxies with a 
substantial young stellar population (H/3> 1.8 A) blur the 
relation and are found at systematically lower Mg b values, 
although over a wide range of FUV— V colour. 

Given the above, is our (FUV— V)-Mgb relation merely 
similar to the original Burstein relation or is it identi- 
ca l? In Figu re [2j we have strived to transform the data 
of BBBFL88 to our own parameters. First, using the IUE 
spectra of BBBFL88, we calculated the 1550— FUV colour 
of NGC 221 and NGC 4649. To c alculate the GALEX 
FUV magnitude of every galaxy in BBBFL88's sample, 
we then simply assumed that each galaxy had one of two 
spectral shapes: either a normal early-type galaxy spec- 
trum like NGC 221, or a UV upturn spectrum like that of 
NGC 4649. The GALEX FUV magnitudes were then sim- 
ply obtained by FUV=1550-(1550-FUV) NGC221.4649, Us- 
ing (1550-FUV)ngc22i for galaxies with (1550- V) > 3.4 
(normal galaxies) and (1550— FUV)ngc4649 for galaxies with 
(1550-I Q < 3.4 (U V upturn galaxies). As the IUE aper- 
tures of IBBBFL881 are not that different from our own, we 
do not correct th e photomet ry for aperture mismatch. Sec- 
ond, we converted BBBF L88|'s Mg? index measureme nts into 
Mg b using the iThomas. Maraston fc Bender! (|2003h stellar 
population synthesis models, assuming an a-element (over- 
abundance [a/Fe]=0.3 for all galaxies. As those Mg2 values 
are central, however, we also need to correct the linestrengths 
for aperture mis match. To do this, we simply compared 
the transformed BBB FL88I Mgfe values with ours for the 
6 galaxies in common (NGC 2768, NGC 4278, NGC 4374, 
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Figure 1. UV upturn dependence on stellar absorption 
linestrengths. The Mg6 (top), Fe5015 (middle) and H/3 (bottom) 
linestrengths are shown as a function of the FUV— V colour, all 
integrated within i? e /2. The top plot is our own version of the 
Burstein relation. Objects inconsistent with a purely old stellar 
population (H/3> 1.8 A) are plotted as crosses, others as filled 
circles. The location of data points lying outside the figure bound- 
aries is indicated by arrows. Note the inverted H/3 axis. Linear fits 
to the data with H/3 ^ 1.8 A are shown as solid lines to guide the 
eye, and the correlation coefficient is reported in the bottom right 
corner of each panel. The best-fit parameters are listed in Tablel2l 
The colour table traces the effective stellar velocity dispersion cr c , 
from purple at small dispersions through green and red at high 
dispersions. 
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Figure 2. (FUV— V)-Mgfo (i.e. Burstein) relation fo r both our 
data (filled circles) and the 1550 A and Mg 2 data of lBBBFL8sl 
(open squares). The latter have been transformed to FUV and 
Mgfe as discussed in the text. Galaxies in common between the 
two samples are connected by a solid line. Error bars have been 
omitted for clarity. The two datasets are consistent, with similar 
scatter. 



NGC 4486, NGC 4552, NGC 5846), and then applied an 
average offset to all th e galaxies ( AMg b= 0.52 A). 

The conversion of lBBBFL88l 's data to our own param- 
eters is necessarily very approximate, but Figure [2] shows 
that our results are consistent with those of BBBFL88. Us- 
ing much improved data and measuring methods, we have 
thus recovered the original Burstein relation. Somewhat dis- 
appointingly, however, the scatter of our data is about the 
same. The scatter is also larger than our measurement er- 
rors, which suggests that it is real. We explore possible rea- 
sons why other studies did not recover the Burstein relation 
in Sections [4] and [5] 



3.2 FUV-V vs. Fe5015 

The middle panel of Figure [1] shows the relation of FUV— V 
with Fe5015. The scatter is very large and a correlation 
analogous to that of FUV— V vs. Mgb is not obviously 
present. This is perhaps not surprising, since the correlation 
of Mgb with the stellar velocity dispersion a is also much 
stronger than that of Fe5015. Equally unsurprising then, to 
our knowledge this correlation has not been discussed much 
if at all in the literature. Interestingly, however, galaxies 
with a young population now appear both below and slightly 
above the best fit to the purely old objects (especially for the 
smaller R e /8 apertures that are not shown here). 

3.3 FUV-V vs. H/3 

The bottom panel of Figure \T\ shows the relation of FUV— V 
with H/3, which is primarily influenced by age. Separating 
young and old galaxies in this diagram is somewhat artificial, 
since we use H/3 itself as our age tracer, but it is striking that 
a tight correlation is again only present amongst the galaxies 
with an old stellar population. 



Having said that, the observed variation of H/3 with 
FUV— V is expected from the metallicity dependence of the 
latter (as traced e.g. by the (FUV— V)-Mg b correlation). It 
thus does not automatically imply an age dependence. 



3.4 NUV-U 

Similarly to Figure [TJ Figure [3] shows the correlations be- 
tween the NUV-U colour and the Mgb, Fe5015 and H/3 
linestrength indices. The contrast between Figure [T] and [3] 
is however striking. First, the NUV— V colour range is only 
about half that in FUV— V for old galaxies. Second, even 
when considering only old galaxies, the UV-linestrength cor- 
relations have essentially disappeared. Third, the galaxies 
with a young stellar component are much more discrepant. 
We also point out that nearly all galaxies with very blue 
NUV— V colours also have high H/3 linestrengths (as ex- 
pected from star formation), although the opposite is not 
always true. 

Although present across the NUV filter, the UV up- 
turn phenomenon usually clearly dominates the spectral en- 
ergy distribution only at wavelengths blueward of 2000 A. It 
thus cannot by itself entirely explain the NUV propertie s ob- 
served. As suggested for exam ple bv lDonas et all (|2007l ; but 
see also iDorman et al][l9 95). who explored slightly differ- 
ent NUV— V correlations, we believe that the above trends 
result from two main competing effects. On the one hand, 
young stars have a greater effect in the NUV than the FUV, 
as only very young stars (t < 0.1 Gyr) can make contribu- 
tions to the FUV light, while young stars in a much large age 
range (t < 1 Gyr) contribute to the NUV. Young stars can 
also influence H/3 for even longer periods (t > 1 Gyr). As tes- 
tified by the (NUV-V")-hU3 plot, this is without doubt the 
reason why the galaxies with a young stellar population are 
so discrepant. On the other hand, line blanketing (the sup- 
pression of the continuum due to a large number of partially 
overlapping absorption lines) is substantial in the NUV and 
increases with metallicity. This is probably why the FUV 
correlations largely disappear in the NUV even for galaxies 
with no sign of significant recent star formation, as early- 
type galaxies generally have a high metallicity. Those com- 
bined effects probably also explain the lack of correlations in 
studi es which lack a true FUV filter (e.g. iDeharveng et all 
120021 ) . 



3.5 FUV-NUV 

Analogous to Figures [T] and [3l Figure [4] shows the correla- 
tions between the FUV-NUV colour and the Mgb, Fe5015 
and H/3 line strength indi c es. In terestingly, as pointed out 
especially by iDonas et all (|2007l ), the UV-linestrength cor- 
relations are actually much more pronouced and show less 
scatter when the UV— optical colour FUV— V is replaced by 
the pure UV colour FUV-NUV. As the metallicity increases, 
this is most likely due to an increase in the FUV upturn con- 
spiring with increasing line blanketing in the NUV. However, 
given that the NUV atmospheric opacity is dominated by 
iron-peak elements, it is surprising that the correlation with 
Fe5015 is poorer than that with Mgb. As expected, how- 
ever, the FUV-NUV colours of galaxies with a significant 
young stellar population are largely uncorrelated with the 
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Figure 3. Same as Figure [T] but for the optical— UV colour 
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0.0 0.5 1.0 1.5 2.0 

FUV-NUV (mag) 



2.5 



Figure 4. Same as Figure [T] but for the pure UV colour 
(FUV-NUV). 



linestrengths. The parameters of the best-fit linear relations 
are again listed in Table [2] 

We stress here that the (FUV-NUV)-Mgfe correlation 
is tighter than the Burstein relation, which uses FUV— V. 
Although the Burstein relation has traditionally been used 
in most work , presumably for historical reasons dating back 
to BBBFL88. and it remains important, our results suggest 



that (FUV-NUV)-Mgb should instead be the correlation 
of choice. As FUV and NUV are generally measured with 
the same instrument, using FUV-NUV rather than FUV— V 
also avoids potential uncertainties in the UV to optical flux 
calibration. 
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Figure 5. Same as FigureQ]but for the FUV-V, NUV-V and 
FUV— NUV colours as a function of the effective stellar velocity 
dispersion <r e . 



3.6 Stellar velocity dispersion 

The Burstein relation and other correlations highlighted in 
this paper, combined with the well-known Mg-cr relation 



fe.g. rTerlevich et al.lll98ll; iBender. Burstein fc Faberlll993l ; 
IColless et al.ll 1999I T Bernardi et al.ll2003f ). imply a correlation 
of the UV-upturn phenomenon with the stellar velocity dis- 



persion, which is a good measure of a galaxy's gravitational 
potential depth and total mass. The colour scheme in Fig- 
ures [TJ [3] and [4] allows us to gauge this dependence, but only 
very roughly. It is thus shown explicitly in Figure [S] The cor- 
responding figure showing the dependence of the UV-upturn 
on the dynamical mass M dyn = 5.0 G _1 Re ol (see lPaper IV|) 
is essentially identical. This mass represents Md yn ~ 2Mi«, 
where Mr/ 2 is the total mass within a sphere containing half 
of the galaxy light. 

Comparing Figures [T] and [S] one can see that the 
scatter in the Burstein relation is smaller than that in 
the (FUV— V)—a e relation. This effect first pointed out by 
BBBFL88 lead them to suggest that the UV upturn is pri- 
marily driven by stellar population properties rather than 
structural or dynamical ones. We confirm this trend here. 
The effect is however much more pronounced when consider- 
ing the FUV— NUV colour, where the correlation with Mgfe 
is very tight but that with a c is much looser. Just as for 
Mg&, NUV— V does not show any clear depend ence on cr e . 

As expected from cosmic downsizing (e.g. ICowie et al.l 
ll996l ; lTreu et alJl2005l ; lThomas et alj|2005h . all galaxies with 
evidence of recent star formation (high H/3) are in the 
low mass half of our sample (<j c < 180 km s" 1 ; see also 
iPaper XVh . 



3.7 Aperture effects 

Analogous UV-linestrength correlations to those discussed 
in this section are also present when using matching elliptical 
apertures of Re/8 and R e . The linestrengths and broadband 
colours are simply slightly offset, as expected from the radial 
stellar populations gradients present in early-type galaxies. 
However, some correlations are clearer for smaller apertures. 
This is most likely because stellar age gradient s are gen- 
erally most pronounce d in gala xy centres (e.g. Paper Vl 
iMcDermid et all l200rj . hereafter IPaper Vllj IPaper XVm . 
The integrated UV— V colours and linestrengths thus have 
much broader ranges of values for small apertures. 



3.8 Multi-dimensional relations 

The various correlations shown in Figures [l] [5] suggest that 
Mg& is the dominant parameter affecting the UV flux, but 
we now consider relations involving three parameters, anal- 
ogously to the Fundamental Plane for stellar dynamics (re- 
lating the stellar luminosity, surf ace brightness and centra l 
velocity dispersion o f galaxies; see lDiorgovski fc Davis|[l987l ; 
iDressler et ailll987h . Without an a priori physical founda- 
tion, such as the virial theorem for the Fundamental Plane, 
it is however not clear what correlations one should expect. 
To minimise biases arising from the very different definitions 
and thus dynamic ranges of the parameters considered (e.g. 
colours, velocity dispersion and linestrengths), we normalize 
each parameter to yield an observed range of 1.0. For exam- 
ple, the measured values of FUV— V extend from roughly 5.4 
to 7.2 mag (yielding an observed range of 1.8 mag), so in the 
fits below N(FUV— V)— 0.0 mag is actually 1.8 mag bluer 
(UV stronger) than N(FUV— V)— 1.0 mag, where the pre- 
fix N denotes a normalised quantity. Likewise for the other 
quantities (FUV-NUV= 0.4-1.8 mag, log(cr c /km s -1 ) = 
1.9-2.4, Mgb= 3.4-4.7 A, Fe5015 = 3.8-5.1 A). 
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Figure 6. Edge-on views of the best-fit planes for selected three-prameter relations (see text). Top-left: Plane relating FUV— V, Mgfe 
and a e . Bottom-left: Plane relating FUV-^, Mgfe and Fe5015. Top-right: Plane relating FUV-NUV, Mgfe and cr c . Bottom-right: Plane 
relating FUV— NUV, Mgfe and Fe5015. The solid lines are shown to guide the eye and all have slope unity and go through the origin. 
The correlation coefficient is reported in the bottom-right corner of each panel. The three-parameter correlations shown are not tighter 
than the two-parameter correlations shown in Figures [T] and [4] 



We have explored a number of parameter combinations, 
but we present only the most promising ones below and in 
Figure [6] in the sense that the resulting planes (in their 
respective three-dimensional spaces) are the tightest. The 
planes were obtained by minimising the square of the resid- 
uals in the direction perpendicular the planes, but we warn 
that the coefficients obtained do depend sensitively on the 
exact quantity that is minimised. The edge-on projections 
of the planes are shown in Figure [6] By construction, the 
solid lines shown (of slope unity and going through the ori- 
gin) should go through the data points (they do). As before, 
the fits are performed to the quiescent galaxies only (i.e. 
H/3sC 1.8 A), yielding 

N(FUV- V) = -3.30 N(Mgb) + 2.00 N(loga c ) + 1.19 , (1) 
N(FUV - V) = -1.95 N(Mg b) + 0.85 N(Fe5015) + 1.21 , (2) 

N(FUV-NUV) = -2.46 N(Mg 6) + 1-25 N(logffe) + 1-07 ,(3) 

N(FUV-NUV) = -1.55N(Mgfe)+0.47N(Fe5015)+1.10 .(4) 

The correlation coefficients reported in the bottom right 
corner of each panel of Figure [6] are no better than those of 
the good two-parameter fits shown in Figures [T] and [3] None 
is better than the (FUV-NUV)-Mgb correlation. This im- 
plies that the FUV excess is not particularly better described 
by employing an additional parameter. This is in agreement 
with our finding of § 13.61 that it is stellar population rather 
than dynamical properties that govern the UV strength. 



4 STAR FORMATION 

In the (FUV-U)-linestrength and (FUV-NUV)- 
linestrength relations discussed in Section [3] the apparent 
correlations defined by the bulk of the observations 
are weakened by the presence of a number of outlying 
data points. However, the outliers systematically have 
H/3 > 1.8 A, which suggests that the deviations from the 
bulk are due to the presence of a young stellar component 
in these galaxies. Indeed, while often negligible in the 
optical, even a small amount of recent star formation can 
have a dramatic effect in the UV. We test this hypothesis 
thoroughly in this section. 

Figure [7] shows how the (FUV-NUV)-Mgfc correla- 
tion measured within R a /2 elliptical apertures depends on 
a number of star formation tracers, none of which is how- 
ever perfect. Our best-fit (FUV— NUV)-Mgfe relation from 
Section 13.51 is overplotted as a dashed line for reference. 
The Ha emission line and the mid-infrared (MIR) emission 
from dust and polycyclic aromatic hydrocarbons (PAHs) are 
normally considered good tracers of respectively unobscured 
and obscured ongoing star formation, while the luminosity- 
weighted mean stellar age traces recent and older star for- 
mation events. The molecular gas mass Mh 2 measures the 
fuel reserve for current and future star formation. As we 
must rely on external sources for those quantities, they are 
not measured within our standard R e /2 aperture, but this 
is unimportant for the argument. We only use those quan- 
tities to gauge the likelyhood of young stars; the data plot- 
ted remain measured within identical apertures. The Ha lu- 
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Figure 7. Dependence of the (FUV— NUV)-Mg b correlation on recent and ongoing star formation. The top panel shows the correlation 
for all galaxies using elliptical apertures of R c /2. Lower panels show the correlation as a function of the integrated Ho luminosity, the type 
of 8.0 (im non-stellar emission if present, the luminosity-weighted mean stellar age and the molecular gas mass normalised by the £?-band 
luminosity. From left-to-right, the columns respectively show galaxies with no, mild or ambiguous, and strong evidence for ongoing or 
recent star formation. A clear correlation with no outlier is only observed for quiescent galaxies. Filled symbols indicate galaxies with at 
least one orbit of GALEX integration time, while open symbols indicate the shorter exposure AIS data. Error bars have been omitted for 
clarity and our best-fit (FUV— NUV)-Mg b relation is ovcrplottcd as a dashed line for reference. 
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mino sity is in tegrated within the SAURON FOV and comes 
from iPaper Vl . In the MIR, we use the Spifczer/Infrared Ar- 
ray Camera (IR AC) 8.0 /im non-stellar emission calculated 
and discussed in IPaper XVl . Considering the emission mor- 
phology, the radial gradient of the 8.0-to-3.6 /J,m flux ratio, 
and the ffpifeer/Infra red Spectrograph (IRS) spectrum when 
available. IPaper XVl verv carefully separated the 8.0 /xm non- 
stellar emission originating from star formation and other 
processes. The mean stellar age is derived through single 
stellar population synthesis models from linestrengths inte- 
grated within a circular aperture of Re/8 (|Paper XVlj) . The 
distance-independent Mu 2 /Lb ratios are global values and 
com e from the CO me asurements of ICombes et al.l (2007) 
and lYoung et all (|2010T) combined w ith the RC3 optical cat- 
alogue jde Vaucouleurs et al.l Il99lf ) . 

The main conclusion from Figure [7] is that a clear, 
tight (FUV— NUV)-Mg b correlation without outliers is only 
present for galaxies which show no sign of current or recent 
star formation. The scatter is much increased for galaxies 
with even mild evidence of star formation, and essentially no 
correlation is present for galaxies with strong evidence. This 
is particularly nicely illustrated by the mean luminosity- 
weighted stellar age panels, which show increasing scatter 
with decreasing mean age. 

Other star formation tracers not shown in Figure [7] 
display analogous behaviours, such as the H/3 absorption 
linestrength and NUV— V colour discussed in this paper, the 
B—V colour, and the 8.0 /jm non-stellar luminosity and asso- 
ciated star formation rate. More importantly, similar results 
are obtained for other apertures (Re/8 and R e ), FUV— V 
and NUV— V colours, and Fe5015 and H/3 linestrengths. 
In addition, the (FUV-U)-(FUV-NUV) and (FUV-V)- 
(NUV— V) colour-colour relations clearly show that, while 
the quiescent galaxies (H/3^ 1.8 A) do define clear correla- 
tions (see, e.g., Boselli et al.l l2005: Donas et al. 2007), galax- 
ies with evidence of young stars (H/3 > 1.8 A) systematically 
depart from them. This shows beyond doubt that most de- 
viations from the correlations described in Section [3] are due 
to galaxies with either current or recent star formation, and 
that the correlations themselves are essentially due to quies- 
cent galaxies (i.e. purely old stellar populations). 

One should not be overly concerned that, occasionally, 
galaxies with apparent signs of current or recent star forma- 
tion in Figure [7] do fall within the region expected of our 
best-fit (FUV— NUV)-Mg6 relation. As mentioned above, 
none of the tracers is perfect, es pecially in e arly-type galax- 
ies. In particular, as shown in IPaper XVj . emission from 
Balmer lines such as Ha in early-type galaxies is normal- 
ity dominated by non-star-formation processes, specifically 
hot evolved stellar populations (PAGB stars). Absence of 
Ha emission thus guarantees absence of star formation, but 
the presence of Ha emission does not imply the presence of 
star formation. Dust and PAHs can also be heated/ionised 
by a variety of proc esses, includ ing these same evolved stel- 
lar populations (see IPaper XVl) . A (FUV-NUV)-Mgfo cor- 
relation is only present for galaxies with 8.0 /im non-stellar 
emission if that emission has a non-star formation origin. 

We also note that the scatter in our correlations is 
not caused by the presence of active galactic nuclei (AGN). 
First, our measurements are global rather than nuclear. Sec- 
ond, there are only four strong radio continuum sources 
in the SAURON sample: two FR I sources (NGC 4374 and 



NGC 4486), NGC 4278 and NGC 4552. None is particularly 
discrepant in our correlations. It is clear however that AGN 
could contribute to the scatter in samples of galaxies con- 
structed differently. 



5 SDSS SAMPLE 

5.1 Sample selection and data 

We now place our findings based on observations of SAURON 
early-type galaxies in the contex t of previous wo rk on 
SDSS galaxies, such as that of iRich et all l|2005h . UV 
and linestrength measurements of SDSS early-type galax- 
ies present a number of challenges. First, most galaxies are 
invariably unresolved by GALEX but the spectral informa- 
tion is derived from SDSS fibre spectra with a 3" diameter, 
leading to aperture mismatch. Second, except at the high- 
est luminosities, many early-type galaxies in the low-redshift 
regime (z ~ 0.05-0.1) are either undetected or barely de- 
tected at the depth of the GALEX Medium Ima ging Sur- 
vey (MIS: iMartin et al.|[2005l ; iMorrissev et alj|2005h . Third, 
a reliable morphological classification is far more challeng- 
ing. Fourth, disentangling the UV upturn from residual star 
formation is very difficult. 

Given that we do find a correlation between UV— optical 
colours and Mg b in SAURON galaxies, we nevertheless attempt 
here to overcome the se hurdles a nd show that the Burstein 
relation pioneered by IBBBFL88I can also be found in more 
distant early-type galaxies, as can other UV-linestrength 
relations. We use for this purpose the M orphologically Se- 
lected Ellipticals in SDSS sample (MOSES: [Schawinski et all 
2007), which includes all spectroscopic SDSS galaxies with 
r < 16.8 in the redshift range 0.05 < z < 0.1 and provides 
reliable visual morphologies. 

The SDSS spectra were first processed with the soft- 
ware GA NDALF to remove nebular emission lines (see 
IPaper Vh. while using pPXF to fit the stellar kinematics 
(jCappellari fc Emselleml I2004T ). The cleaned spectra were 
then degraded to the spectral resolution of the Lick indices 
and stellar absorption linestrength s were measured, i nclud - 
ing Mgfo and H/3. For d etails, see Schawins ki et al. | <|2007h 
and lThomas et"aL (2008). With those data in hand, we now 
address the four challenges listed above. 

First, there is no reasonable correction we can apply to 
the issue of aperture mismatch, so this necessarily remains 
a source of systematic error for GALEX-SDSS galaxies. 

Second, since the slope of the Burstein relation is not 
particularly steep, there is a danger that it could be washed 
out by large errors in the UV photometry. Even though these 
requirements substantially reduce the sample size, we thus 
only consider here those MOSES early-type galaxies that are 
detected by GALEX in both NUV and FUV filters, and we 
further require that the photometric error in each be less 
than 0.3 mag. These two requirements reduce the sample 
from 12 828 galaxies for MOSES to 535. Perhaps more impor- 
tantly, these requirements strongly bias our sample in favour 
of UV-blue galaxies, particularly compared to the SAURON 
galaxies discussed in this paper. UV-red galaxies are largely 
excluded. 

Third, it is clear that any automated early-type galaxy 
selection criteria based on structural parameters leads to the 
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inclusion of many disc-dominated galaxies, Sa galaxies be- 
ing a particularly prevalen t source of contamination (see, 
e.g.. lSchawinski et alj|2007h . As galaxies with discs are likely 
forming stars, they will contaminate the Burstein and other 
similar relations with UV flux from young stars. The most 
reliable method for excluding such non-early-type interlopers 
is visual inspection, which has been carefully performed for 
the MOSES sample using SDSS images. Short of obtaining 
space-based data for a substantial number of these galaxies, 
this is thus the best morphological classification possible. 

Fourth, even after removing intermediate and late- 
type interlopers, residual star formation i s still common 
in low-redshift early-type galaxies (e.g . lYi et aL 2005; 
ISchawinski et af]|2007l ; lKavirai et al.ll2007l ). and this residual 
star formation will dilute the Burstein and other relations. 
To remove these objects from the sample, it is possible to re- 
ject all galaxies with high H/3 linestrengths. Specifically, we 
investigate here rejecting objects with H/3 + H/3 Grr > 1.8 A 
(see below). 



5.2 Burstein relation: FUV-V vs. Mgb 

We plot in the top panel of Figure [8] the Mg b linestrength 
against FUV— g colour for the SDSS early-type galaxies se- 
lected. Overplotted is the best-fit Burstein relation derived 
in Section f3.ll transformed from V to g. The latter correc- 
tion was obtained using a large sample of early-type galax- 
ies with GALEX data selected visually from SDSS. As the 
g — V colour effectively does not vary with FUV— NUV, we 
adopted a fixed correction g — V = 0.53 mag. 

Ignoring the cut on H/3 linestrength, the data in the 
top panel of Figure [8] only show a loose correlation in a 
sense opposite to the Burstein relation discussed in Sec- 
tion 13.11 The immense scatter introduced by residual star 
formation largely washes out the signal in the UV— optical 
colour. If we however separate the sample according to H/3 
linestrength, the picture changes dramatically. Quiescent 
early-type galaxies with low H/3 values (161 objects with 
H/3+H/3err ^ 1.8 A; black circles) cluster at redder colours 
and are located roughly within the region expected from the 
original Burstein relation, with few outliers. Galaxies with 
a young population with high H/3 values (144 objects with 
H/3-H/3 e rr > 1-8 A; grey circles) do not and have an entirely 
different distribution. This selection on H/3 does exclude 
galaxies whose error bars bestride the H/3= 1.8 A thresh- 
old between recently star-forming and non-recently star- 
forming galaxies (i.e. galaxies for which H/3-H/3 crr ^ 1.8 A 
< H/3+H/3 crr ), but we believe this is best suited to test our 
hypothesis that star formation (i.e. the presence of young 
stars) is mainly responsible for the outliers in the Burstein 
relation. Galaxies that bridge the threshold may or may not 
contain young stars and will blur any dichotomy. 

Our H/3 criteria are however not perfect discriminants 
between entirely passive early-type galaxies and those with 
a low level of residual star formation in the last gigayear. In 
any case, residu al star formation can occur at large radii 
(see iPaper XHlf) . and this will not be picked up by the 
SDSS fibres. If we use a blunter selection criterion and re- 
strict ourselves further to early-type galaxies with a high 
velocity dispersion (a > 200 kms -1 ; 169 of 535 galaxies), 
where r esidual star formation is virtually completely sup- 
pressed l|Schawinski et all 120061 . see, e.g.,), then we can see 
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Figure 8. Mg-(FUV-g) relation for the MOSES sample early- 
type galaxies with good UV photometry. Top: Entire sample meet- 
ing the H/3 selection criteria. Bottom: Only galaxies with a ve- 
locity dispersion a larger than 200 km s — 1 . Quiescent galaxies 
with H/3+H/3 crr ^ 1.8 A are shown as black circles, while galax- 
ies consistent with harbouring a young stellar population (H/3- 
H/3 err > 1.8 A) are shown in grey. Characteristic error bars are 
shown in the top-left of each plot. Our best-fit Burstein relation 
is overplotted as a dashed line for reference, taking into acount 
the small g — V colour transformation described in the text. 



that all early-type galaxies cluster within the same locus re- 
gardless of the H/? criteria (see Fig. bottom panel). And 
indeed, there are very few galaxies with H/3-H/3 orr > 1.8 A 
left (5 objects compared to 97 for H/3+H^ crr < 1.8 A). The 
large errors in both Mgb and FUV— g, in addition to aper- 
ture mismatch, result in sufficient scatter to largely hide the 
(FUV— g)— Mg b relation, but the locus is consistent with the 
slope and zero-point found for the SAUR0N galaxies. 

Our treatment of the SDSS sample thus supports the 
presence of (UV— optical)-linestrength relations as described 
in Section f3. 11 but again only for quiescent galaxies. This is 
true using both an H/3 linestrength (i.e. age) cut and a stellar 
velocity dispersion (i.e. total mass) cut. Our treatment also 
illustrates how challenging weeding out galaxies with low 
levels of residual star formation can be with the data quality 
typical of large surveys. Here, the sheer number of galaxies 
does not trump the need for high quality data. 
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6 DISCUSSION 



6.1 Origin of the UV upturn 



A main result of our work illustrated in Figures [THI] is that 
we recover the Bu rstein relat ion for early-type galaxies origi- 
nally discussed in lBBBFL88l . as well as the (NUV-U)-Mgfe 
and (FU V— NUV)-Mg6 correlations discussed in later pa- 



pers (e.g. Dormanetal 



_ J[ 19951 ; [Rich et al.ll2005l ; iBoselli et ail 
120051 ; [5 onas et al.l 120071 ). However, contrary to these au- 
thors, we use identical apertures for all quantities. We 
have also clearly showed that correlations with the H/3 
linestrength index also exist, albeit usually with larger 
scatter. Correlations with Fe5015 are at best very weak. 
These additional correlations are crucial, as once compared 
to stellar population synthesis models they will allow us 
to roughly separate the effects of a-element enhancement, 
metallicity a nd age on the UV-upturn phenomenon (see, e.g., 
iPaper XVIl|) . 

Lacking other elements, previous wo r ks considered Mg 
as a metall i city tracer (e.g. IBBBFL88I ; I Rich et al.l 120051 ; 
iDonas et al.lf2007h . But given that early-type galaxies and 
bulges are known to be over-abundant in a elements (e.g . 



Rich! 19881; Gorgas. Efstathiou fc Aragon Salamancal ll99(J ; 



Worthev. Faber fe Gonzalej Il992l )." one might now naively 
be tempted to treat Mg primarily as a tracer of a-elements, 
and take Fe5015 as the prime metallicity tracer. A high 
[a/Fe] ratio is usually interpreted as implying a short star 
formation timescale (i.e. a burst), as stars are enriched in 
a elements produced by early Type II supernovae but are 
not enriched in Fe produced by longer t imescale Type la 
supernovae (e.g. lEdvardsson et all 1 19931 ; iFuhrmannl Tl998; 
Ffh omas. Greggio fe Bended 1 199&t ). Taken at face-value, the 



tight (FUV— V)-Mg b correlation (i.e. the Burstein relation) 
would thus suggest that the UV upturn in early-types is pri- 
marily due to stars formed in a rapid burst and enriched in 
a elements. The extremely loose (FUV— U)-Fe5015 correla- 
tion would further suggest that enrichment in Fe produced 
over longer timescales is not the primary driver of the cor- 
relations. This would then presumably pose challenges to 
UV upturn theoretical models with strong metallicity de- 
pendences (e.g. EHB models), and equally spell bad news 
for s cenarios relying on binary star evolution m ore generally 
(e.g. iMaxted et aT1l200ll ; lHan et alJkool 120071 ) . 

However, are we really justified to make this simplistic 
association of Mg6 with a elements and Fe5015 with metal- 
licity? The short answer is (un)fortunately no. Although Mg 
is an a element, the Mg6 index is still primarily driven 
by metallicity, while the Fe5015 index includes a number 
of poorly understood absorption lines. Furthermore, while 
they slightly affect both, a elements dominate neither the 
energy generation (He, CNO) nor the atmosphe ric opacity 
(iron -peak elements) in metal-rich galaxies (e.g. iGuo et all 
l2010h . Beside specific spectral featu res, thei r effect on UV 
light should thus be limited (see also lY"l[2008l ). A proper un- 
raveling of the individual effects of a-element enhancement 
and metallicity on the UV upturn is thus carried out in the 
companion paper by Jeong et al. (in preparation), with the 
help of the s tellar population synthesis models presented in 
IPaper XVlfl , This is not only the case for the primary de- 
pendence of the UV upturn on metallicity, but also for its 
dependence on H/3 for old galaxies (H/3 ^ 1.8 A). A priori, 
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Figure 9. Same as Figure [T] but for the FUV— NUV colour as a 
function of the specific angular momentum of the galaxies . The 
normalised number distribution of slow-rotator (red) and fast- 
rotator (blue) galaxies as a function of the FUV— NUV colour 
is shown on the right for quiescent galaxies only (i.e. cirles and 
triangles). 



this may well be explained entirely by metallicity effects, but 
it could also hide a residual age dependence. 

Despite having failed to find a dominant dependence of 
the UV upturn on <r c , we can test if other parameters might 
better reveal a structural or dynamical influence on the UV 
upturn. We find no or weak dependences of the FUV— V, 
NUV-V and FUV-NUV colours on the ellipticity e, boxi- 
ness of the isophotes 0,4 and the ratio of rotational to pres- 
sure support V/a. In the spirit of recent results of the SAUR0N 
survey, we plot in Figure [§] the dependence of the UV up- 
turn on Ajj c , the speci fic angular momentum of a galaxy 
measured within R c (see lPaper IXl for details). Figure [9] also 
shows the normalised number distribution of slow-rotator 
(A_r c < 0.1) and fast-rotator (Xr c 0.1) galaxies as a func- 
tion of the FUV— NUV colour, for quiescent galaxies only. 
Although there is much scatter, we see that among quiescent 
galaxies very blue FUV-NUV colours are more commonly 
observed among slow rotators than fast rotators. Clearly, 
however, more extensive investigations are necessary to sub- 
stantiate any dynamical dependence of the UV upturn. 

Lastly, we point out that the scatter in the UV- 
linestrength relations presented here (e.g. the Burstein re- 
lation) is larger than our measurement errors. It is thus real, 
implying that there is not a one-to-one correlation between 
the UV excess and absorption linestrengths. This is an im- 
portant result, but stellar evolution models of course do not 
depend on linestrengths but rather on age, metallicity and 
a-element abundance. 



6.2 Star formation 

Our work shows that star formation is responsible for most 
outliers in the UV-linestrength correlations. Indeed, with 
the extensive ancillary data available for the SAUR0N sample, 
and the exquisite sensitivity reached for nearby objects, we 
were able to systematically identify all outliers with galax- 
ies showing signs of recent or current star formation. This 
could also be done in an analogous manner for a much larger 
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(and more distant) sample of galaxies from SDSS. The fact 
that clear correlations are then only recovered for massive 
galaxies with high quality data can probably be explained 
by the relative poor sensitivity of SDSS to the presence of 
recent and/or current star formation (more likely present in 
low-mass systems). 

The s a me r eason probably explains the results of 
Rich et al.l (|2005r ). which are most at odds with those of 
BBB FL88I . We support here the results of iDonas et al.l 
(|2007l ), who argued in favour of contamination from low- 
level star formation in their own sample, especially for the 
lenticulars. Requiring S/N^ 10 for the SDSS spectra (al- 
though t hey quote a final S/N= 23 ± 7 for their quiescent 
sample) . iRich et al.l |2005l ) were probably not as sensitive to 
weak emission lines and thus star formation as we are here. 
It would be interesting to see the UV— optical colour cor- 
responding to their Ha emission (and thus star formation) 
upper limits. We require here S/N^ 60 per spatial bin (and 
spectral element) for SAURON and simultaneously fit both ab- 
sorption and emissio n lines, yi elding an equivalent width sen- 
sitivity of 0.1 A (see lPaper Vh . We also have much ancillary 
data, including the age-sensitive H/3 absorption line index. 

We also stress here a fact that has generally been over- 
looked in previous UV-upturn observational studies. All 
galaxies which are off the relations, and which we identify as 
having current or recent star formation in Sections U and 
are systematically below the correlations. That is, they sys- 
tematically have Mgb and Fe5015 linestrengths which are 
lower than those predicted by the correlations for purely old 
galaxies. None has significantly higher Mg6 or Fe5015, and 
the effect is particularly pronounced for Mg6. Considering 
the predictions of s tellar population synthesis models (e.g. 

omas et al.ll2003T ). this is entirely consistent with the out- 
liers being of systematically lower mean ages. 



7 CONCLUSIONS 

We have used space-based UV ultraviolet photometry from 
GALEX, ground-based optical photometry from MDM and 
ground-based optical integral-field spectroscopy from SAURON 
to study the UV-linestrength relations of early-type galax- 
ies. Focusing on the SAURON sample, those data represent 
significant improvements over previous work. In particular, 
identical apertures could be used for all datasets, eliminat- 
ing the aperture mismatch issues that have affected previous 
studies. 

We have shown that early-type galaxies follow fairly 
tight correlations between the integrated FUV— V and 
FUV-NUV colours and the integrated Mg& and H/3 
linestrength indices, less so for the NUV— V colour, Fe5015 
index and stellar velocity dispersion a. Through stellar pop- 
ulation synthesis models, these correlations constrain the 
dependence of the UV-upturn phenomenon in early-type 
galaxies on age, metallicity and Q-element abundance. In 
particular, despite some recent controversy stemming from 
the analysis of large survey datasets, we have recovered a 
(FUV— y)-Mg b correlation that is an alogous an d entirely 
consistent with that first proposed by lBBBFL8a . We refer 
to this and the original relation as the Burstein relation. 
This relation clearly suggests a positive dependence of the 
UV-upturn on metallicity and is stronger when the pure UV 



colour (FUV— NUV) is used, although any possible depen- 
dence on age must await a better comparison with models. 
The scatter in the correlations appears to be real and there 
is mild evidence that a strong UV excess is preferentially 
present in slow-rotating galaxies, although this last point 
deserves further study. 

Using ancillary data constraining past and current star 
formation in our sample galaxies, we have shown that most 
data that do not follow the main correlations can be at- 
tributed to galaxies with evidence for recent star formation, 
even if weak. The UV-linestrength correlations discussed 
thus appear to hold exclusively for old stellar populations. 
A large sample of more distant early-type galaxies selected 
from SDSS, and thus easier to relate to current large surveys, 
supports this view. The same can be said of the correlations 
outliers themselves, which have systematically lower Mg6, 
Fe5015 and higher H/3 linestrength indices. 

Given that galaxies show internal colour and 
linestrength gradients, it is natural to ask whether the 
global correlations observed here also apply locally within 
galaxies. Using the full capabilities accorded to us by 
integral-field data, this issue is the focus of a companion 
paper (Jeong et al., in preparation). 
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APPENDIX A: AVERAGE LINESTRENGTH 
MEASUREMENTS 

Formally, to measure the average linestrength within a 
SAUR0N aperture containing multiple bins, we should sum 
the (fully calibrated) spectra of all the bins within that 
aperture and re-measure the linestrength. However, recal- 
culating the linestrength implies re-deriving the stellar kine- 
matics, re-subtracting the ionised gas emission, re-correcting 
the linest rength for velocity broadening, etc, as done in 
iPaper Vj . It is thus a time-consuming task. Furthermore, in 
a companion paper exploring the UV-linestrength relations 
as a function of radius within individual galaxies (Jeong et 
al., in preparation), we integrate the linestrengths over el- 
liptical annuli (rather than elliptical apertures). The result- 
ing LOSVDs can then be significantly non-Gaussian, which 
makes the stellar kinematic and velocity dispersion correc- 
tion calculations problematic. It is thus preferable to find 
a method to integrate (i.e. average) the linestrengths in an 
aperture using the linestrength measurements in individual 
bins themselves. 

In Figure I All we s how the pr operly averaged 
linestrength measurements of IPaper XVIll within a circular 
aperture of radius _R C , compared to measurements obtained 
by simply averaging the linestrength values of the individual 
bins within that same aperture, weighted by the luminosity 
of each bin. While there is some scatter, it is much smaller 
than the range of values covered, the difference between the 
two measurements is alwa ys small, typi cally less than the er- 
ror on each measurement l|Paper XVlJ quotes a formal ran- 
dom uncertainty on integrated linestrengths of at most 0.1 A, 
and we have adopted the same uncertainty on our measure- 
ments here), and the few outliers systematically have very 
large effective radii. Most importantly, there is no significant 
trend with the linestrength values themselves. An aperture 
of R e is also the worst case scenario, and even smaller differ- 
ences are obtained for smaller apertures involving fewer bins. 
It thus seems that, although not formally equivalent, sim- 
ple luminosity-weighted linestrength averages do yield accu- 
rate estimates. We have thus adopted this averaging method 
for SAURQN papers req uiring int egrated line strength values 
(e.g. JScott et al.ll2009l . hereafter IPaper XlVt Jeong et al., in 
preparation) . 
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Figure Al. Comparison of integrated linestrength values ob- 
tained by properly summing the spectra (labeled "sum" ) and sim- 
ply luminosity-weighting the linestrength values themselves (la- 
beled "Iwa" ) within circular apertures of radii R c . Error bars are 
omitted for clarity, but the adopted formal random error on either 
quantity is 0.1 A. Top: Mgb. Middle: Fe5015. Bottom,: H/9. The 
scatter is always much smaller than the range of values covered 
and there is no significant systematic trend. The dotted lines show 
the 0.14 A formal random error on the difference. 
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